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Dissipation of tidal energy in the moon was calculated under the 
< ssumption that it can be represented as due to imperfect elasticity. If 
the factor l/Q for dissipation per cycle is assumed to be 1/100 for 
distortions strain energy and 1/1000 for dilational strain energy, the 
heat now being generated in the moon by tides is less than .010 ergs/gm/w- = 
i.e., negligible compared to radioactive heating by a chondntic composition. 
Tidal heating would be comparable to radioactive heating, however, if the 
semi-major axis of the moon's orbit was one-third as great, so the mechanism 
does limit the possible history of the moon's orbit. Also, appreciable tidal 
dissipation would cause thermal stresses and be conducive to convection, 
since it is much greater in the center than near the surface of the moon, and 
has a non-uniform distribution in latitude (and longitude as well, if • 

rotation is synchronous wit h revolution). 

Urey a nd others [19593 suggested that tidal dissipation may have 
significantly contributed to heating of the moon's interior in the past. 

Kopal [1962] has calculated the tidal dissipation in a moon assumed to be 
t fluid with Newtonian viscosity. However, since tidal distortion is 
relatively high frequency phenomenon, it seems a more modest extrapolate a 
from experience to assume that tidal dissipation in the moon occurs as, a 

result of imperfect elasticity with a factor l/« for dissipation per cycre 

+-,,3 f 0r -hhe earth's mantle from polar tides, fiee 
comparable to those estimated tor mne etw. uu a 

oscillations, and latitude variation. — . , 
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The tidal disturbing function W is: 

JL 

Where P. is a Legendre Polynomial; Y, X are radius, latitude, and 

V 

longitude in a moon-fixed coordinate system; the asterisked quantities 
refer to the disturbing body; C-M* is tho product of the gravitational 
constant and the mass; and S is the arc from (0*A*) to (0,X). Apply the 



Where is the Legendre Associated Function and 6 Qm is the Kronecker delta. 

To obtain the variation in time of V*,0*A*> they must be expressea 
in terms of the earth's orbit referred to the moon. For the present orbit 
this is most simply done using the numerical values of coefficients as given 
by Brown's theory of lunar motion, as has been recently described by 
Harrison L1963]. Because we wish to investigate the effects of changing the 
orbital parameters, and because the necessary computer subroutines had already 
been written for another purpose (to degree of detail superfluous to the 
present problem), we used a different development [Kaula, 1961] based on the 
assumption that the orbit can be considered a Keplerian ellipse at a 
fixed inclination to the moon's equator with secularly moving node and perigee 
For the tidal problem, the most significant omissions under this assumption 
are short-period perturbations of the semi-major axis and the longitude by 
, he sun. For the present orbit, the largest of these terms (those containing 
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h in the arguments of equations (6)-(8) of Harrison [ 1963 ]) have a ratio 
about 0.2 to the terms arising from ellipticity of the orbit. This ratio 
vould vary directly with variation in the semi-major axis of the lunar orbit, 
but would stay about the same with variation in the inclination or 
eccentricity. So it is consistent with the unavoidable crudeness of our 
estimates of the dissipation factors l/Q to assume a purely elliptic orbit. 
Apply the transformation of equations (7)-(28) in Kaula [ 1961 ] to 
(sin Ceos mX*, sin mX*3 equation (2): 
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Where a*, e*, i*, 0*, u>*, M* are the Keplerian elements of the earth's orbit 
referred to the moon's equator and a departure point thereon fixed with 
respect to inertial space; F^U*) and O^e*) are polynomials of the sine 
and cosine of the inclination and of the eccentricity, respectively, ana 6 
is the "lunar sidereal time" : the angle between the inertially fixed departure 

point and the point on the moon from which selenographic longitudes are 
measured. We assume the rate 0 to be constant, which is equivalent to 


neglecting the physic si lihrs,tion» 
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We abbreviate equation (3) as • 
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Where the single subscript g replaces the subscript pair p,q. 


the amplitude 
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The strain energy per unit volume, divided into shear (or distortional) 
and compressive (or dilatational ) terms [p.12, Jeffrey s, 1959]: 
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summing over repeated subscripts. In equation (\) , , p, is the rigidity, k 
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is the bulk modulus , e ±j is the strain tensor and e'^ is the part of the 
strain tensor expressing departures from symmetry: 
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where is the Kronecker delta. 

If the tidal disturbing function is expressed as a sum of spherical 


harmonics : 




( 9 ) 


the strain tensor e^ at any point (V,0A) can be expressed as 



where €. 1J(Uac ,s is the 
modulus } bulk modulus } 
disturbing function. 


response of a planetary model of specified shear 
and density to a unit coefficient in the tidal 
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Comparing equations (4) and (9), we have: 



Substituting equations (10), (ll), (12) in equation (7) and sorting 
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For any particular term of subscripts 't,,m>u,v,g,h, the energy 
dissipated in one cycle of duration 2n /( o^ + cr uvh ) will be: 



and in one cycle of duration 2 ti/ ( " ^uvh^ * 



(16) 
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To obtain the contributions to enorgy dissipation per unit time, 
multiply equations ( 15 ) and- (l6) by th< absolute values of the rates per 
unit time, | 0<Bg + d u J and | 0<Bg - c uvh I- Sum over all subscripts to obtain 
the mean dissipation rate at a particular point fixed in x,he moon: 
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If the rotation is synchronous with revolution, as exists at present. 
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there are terns which will contribute only through their amplitude a^ g but 
not through their rate a ^ In the synchronous case, a reference longitude 
±ust be fixed. Making this reference longitude the mean direction of the 
earth sets zero all terms containing sin {m(<D* + M* + O* - 6)} and sets unity 
all terms containing cos {m(w* + M* + ft* - 6)}: i.e., in equation (10) the 
contribution to (t) will be a^ g for l-m even and 0 for l-m odd, and to 
%m (t) it will be 0 for <t-m even and for £-m odd. Practically the only 

term for which this effect is significant is (£,m,p,q) = (2,2,0, 0). 
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In addition, there will be degeneracies for m = 0 in all cases and for 
m ^ 0 in the synchronous case, requiring the combination of terms before 
proceeding as in equations (l3)-(l6). In these cases, the rate for the tern 
of subscripts U,m,p,q) will be the negative of the rate for the term of 
subscript U,m, i-m-p, -q). If the rate is taken as that of the U,m,p,q) 
term, then the amplitude for the cosine coefficient will be: 



where the subscript correspondence is h with (p,q) and i with (t-m-p, -<l)> 
and a^, are computed by equation ( 5 ). Then equation (l4) must be 

modified so that b , b^ coefficients appear in front of ^ijuvs' 


etc. terms in place of a^, a uvh . 


Another set of degeneracies occuring 

tctrns of Subset t p t J nrr> } Pi 4-0 

to terms of subscript (£, nH-2i, p-i> 


in the synchronous case arises 

v,{ll bftve Vfttes 

q), where i is any integer. 


Including terns for which the disturbing function rate is zero in 


effect makes the energy dissipation rate a 


function of the constant value of 


the strain, which raises the question of whether strains from other thar tidal 
causes should be considered as well. Since we are interested in the 
dissipation over geological durations of time, these terns perhaps should 
be omitted because in such time we would expect non-oscillating strains to 
he removed by anelastic processes. But then the dissipation rate obtained 
would be an absolute minimum for the assumed Q. The moments of inertia 
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of the moon indicate that it now contains strains larger than tidal, so 
leaving in the non-oscillating tidal terms should yield a dissipation rate 

unlikely to be too high. 

The quantities a^, O ^ l - 2, m = 0,1,2 where calculated from 

equations ( 5 ) and (6), for a variety of lunar orbits, using subroutines 

•For f (i*) and G, (e*) originally devised for analysis of close 
<lmp v ' tpq 

satellite orbits, and computing the rates to*, M*, Ly the methods 
described in Kaula [196l]* 

The strain tensors were calculated using the formulation 

of the earth tide problem of Alteman and others [1959], which has also 
been used by Takeuchi and others [1962] and Longman [1963]. In this 
formulation, the basic variables are the radial factors of vector spherical 
harmonic expressions of the displacements, stresses, and potential terms. 
y of the radial displacement; y 2 , of the compressive stress; v 3 > of the 
tangential displacement; y^, of the shex- stress; y^, of the potential, 
y 6 , of the potential gradient less the contribution thereto of the radial 
displacement. The equations of equilibrium then become a system of six first 

order equations: 



( 21 ) 


The P are functions of V, k, y,, g, and p. Three columns of P.. have terms 
of 0(V 2 ) , so the requirement of regularity at the origin eliminates three 
constants of integration. The three surface conditions of zero tangential 
stress, zero radial stress, and the potential gradient being related to the 
poetntial as a spherical harmonic in free space make the problem determinate 
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After solving equation (21 ) numerically for the y^'s corresponding to a 
particular harmonic = the contribution to the strain matrix 

£ is calculated by (here 9 is colatitude, 0 is longitude, and X is 
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Equation (J&) can be derived using the equations on page 56 of Love [ 1927 ], 
applying a factor of 1/2 to the off diagonal components to be consistent 
with tensor convention. The 6^’s from the Y^'s were used in equation (ik). 
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The pole tide suggests a Q of more than 100 for the earth at a 14 -month 
oeriod [p. lo2 , Munk & MacDonald , I960], while the latitude variations 
Indicate a Q of about 40 at the same period [p.l48, Munk & Ma cDonald, i960; 
pp. 255-259, Jeffreys , 1959]* The free oscillations of the earth indicate Q's 
of 170 to 400 for periods less than an hour r Benioff and others , 196lj 
Ness and others , 1961], except for a Q of 7500 for the S Q °, the only mode 
which is purely compressive [ Ness and others , 1961]* Models for rock creep 
proposed by Jeffreys and Cramp in [i960] «ad J. R. MacDonald [1961] ^suggest 
Q's between 60 and 100 for semi-monthly and monthly periods. A shear Q s 
of 100 thus seems a reasonable compromise. Considering that excitation from 
the atmosphere and oceans may maintain the S 0 ° free oscillation [ Ness and oth e r s, 
1961], we assume a compressive Qc of 1000. 

The strain energies were calculated for several lunar models proposed 
by Harrison [1963]. However, since the uncertainty in Q reduces this problem 
to one of estimating order of magnitude, this discussion will be limited to 
a homogeneous moon of density 3*3^ " m ~3, rigidity 7- 3^ x TO- 1 - 1 dynes cm 2, 

and bulk mod ul us 1.23 x 10^ dynes cm The Love numbers obtained for ..his 
model by the numerical solution of equation (15) were .0344 for h, .0195 f° r 
and .0095 for l. 

The thermal histories of lunar models with chondritic composition have 
been calculated by G. J. F. MacDonald [1959]* Even with a cold origin, these 
models come very close to melting at depths in the moon exceeding 500 xm. 

The chondritic composition used had radioactive contents of 8.0 x 10 , 

1.1 x 10" 8 , and 4.4 x 10~ 8 g/g for potassium, uranium, and thorium, respectively, 
which yield a thermal energy output of 1*59 ergs/g/year at present and 
12.8 ergs/g/year 4. 5 x 10^ years ago. Hence for tidal dissipation to be 
significant, it should contribute on the order of 5 ergs/g/year, or .04 ergs/ 
cm3/ day. 
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The results for the present orbit of the moon are shown in Figure 1 
in the form of maps of the energy dissipation at four levels within the 
moon. Since the dissipation is symmetric about the equator and two 
meridians at right angles, we need show only one octant for each level. 

The evident features are firstly that the tidal dissipation is at present a 
negligible source of heat, and secondly that the distribution of the 
heating is extremely non-uniform both radially and laterally. This 
variability of distribution suggests that if the moon's orbit had ever 
been such that the overall heating was appreciable, large thermal stresses 
would have resulted, leading to convection or some other form of mass motion. 

The amplitudes and rates were therefore calculated for different 
orbital specifications. In turn each element was varied, holding the 
others fixed at the present values: l) the semi -major axis (holding the 

rotation synchronous with revolution about the earth); 2 ) the eccentricity; 

3 ) the inclination; and 4) the rotation rate. The results are displayed in 
Figure 2, in the form of curves showing variation in the average ratio to the 
present dissipation with variation in the orbital elements. The variation 
is particularly marked with variation in the semi-major axis. If the semi- 
major axis were only one-third as great as it is now, the criterion of 
.Ok ergs/cm 3 /day would be exceeded for most of the moon. At the secular 
acceleration calculated by Munk & MacDonald [i 960 ], the moon would have 

been at this distance about 10 years ago. 

The pattern of energy dissipation shown in Figure 1 is composed of 

even degree harmonics symmetric about the equator: (4,4), (4,2), (4,0), 

(2,2), and (2,0). In a moon close enough that heating by tidal dissipation 
was large enough to cause convection, the 2 d degree terms in its mass 
distribution would be determined mainly by the gravitational attraction of 
the earth. However, the fourth degree terms would be determined 
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by the convective pattern. It therefore will he of interest when 
variations in the gravitational field and external form of the moon are 
better determined, to find out whether these fourth degree terms are 
markedly larger than other terms, such as the third degree terras. If they 
are, it would be a strong indication that the moon was once close enough 
for heating by tidal dissipation to cause convection, and hence would lend 

further evidence as to the moon’s origin# 

In conclusion, we can say that heating by tidal dissipation is currently 
insignificant (unless Q factors from the earth's mantle are wrong by a 
factor of about 10 2 ), but that this study confirms the suggestion of Urej£ 
end others [1959] that it may have been important in the past when the moon 
came much closer to the earth. 
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